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ABSTRACT 

A theory of transport coefficients in weakly compressible turbulence is derived by applying 
Yoshizawa’s two-scale direct interaction approximation to the compressible equations of 
motion linearized about a state of incompressible turbulence. The result is a generalization 
of the eddy viscosity representation of incompressible turbulence. In addition to the usual 
incompressible eddy viscosity, the calculation generates eddy diffusivities for entropy and 
pressure, and an effective bulk viscosity acting on the mean flow. The compressible fluc- 
tuations also generate an effective turbulent mean pressure and corrections to the speed of 
sound. Finally, a prediction unique to Yoshizawa’s two-scale approximation is that terms 
containing gradients of incompressible turbulence quantities also appear in the mean flow 
equations. The form these terms take is described. 
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I. Introduction 


This paper derives a gradient transport model for weakly compressible turbulence 
which generalizes the eddy viscosity description of incompressible turbulence. The model 
is derived by applving Yoshizawa s two-scale direct interaction approximation (TSDIA) 
to the theory of weakly compressible turbulence. In this theory, the compressible field 
quantities are analyzed by linearization about a state of incompressible turbulence. 

The present model can be compared to Yoshizawa’s theory , 4 the first comprehensive 
attempt to derive a general transport model for compressible turbulence. Like Yoshizawa’s, 
the present theory requires the derivation of a direct interaction approximation 5 (DIA) for 
compressible turbulence. Whereas Yoshizawa’s theory allows strong compressibility effects, 
the restriction in the present work to weak compressibility makes possible a somewhat more 
complete treatment of the DIA response functions and leads to more explicit expressions 
for the transport coefficients. In particular, the present theory treats the complete matrix 
of response functions, including various coupling terms which produce the wave motions 
characteristic of compressible flow. 

The DIA for weakly compressible turbulence derived here can be compared to the- 
ories proposed by Bertoglio et al . 6 Like them, we derive approximate formulas for the 
response functions; although a comprehensive DIA theory of compressible turbulence was 
formulated by Hartke et al , 7 explicit expressions for the DIA field descriptors were not 
given. Renormalization group methods were applied to strongly compressible turbulence 
by Staroselski et al . 8 The spectral dynamics of weakly compressible turbulence has also 
been modeled using EDQNM by Bertoglio et al 8 and by Bataille et al. 

The theory of weakly compressible turbulence has been considered from the TSDIA 
viewpoint bv Shimomura ; 10 however, this work had the quite different goal of a more refined 
treatment of buoyancy effects than is possible in the Boussinesq approximation. Another 
refinement of the Boussinesq equations through low Mach number expansions is investi- 
gated by Mlaouah et al . 11 A distinct theory of weak compressibility has been proposed by 
Girimaji . 12 In this theory, the compressible part of the velocity field is decomposed into 
pressure and temperature-dominated components. Any of these theories could have been 
used as a starting point for TSDIA instead of the straightforward linearization adopted 
here. 
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A more traditional approach to modeling compressible turbulence generalizes incom- 
pressible models through mass-weighted averaging; a definitive survey of this viewpoint 
has been given by Huang et al . 13 Like Yoshizawa, we attempt instead to derive a model 
from a dynamic theory of compressible turbulence. Ristorcelli 14 has proposed a pseudo- 
sound theory of compressible turbulence for compact flows of engineering interest. In this 
weakly compressible theory the length scales of the acoustic fluctuations scale as the in- 
verse Mach number with respect to the vortical length scales. This allows the neglect of 
wave propagation effects. 

The present model is derived from a three-component decomposition of the com- 
pressible flow field: after the usual Reynolds decomposition into mean and fluctuating 
components, the fluctuations are further decomposed following Zank and Matthaeus 2 and 
Erlebacher et al 3 into an incompressible field and small compressible perturbations. The 
compressible field is treated by linearization about the incompressible field; it follows that 
this theory only applies to compressible turbulence at small turbulent Mach number. 

As usual, the Reynolds averaging introduces various correlations of fluctuating quanti- 
ties which require closure; whereas for incompressible turbulence, only the Reynolds stress 
arises this way, the increased number of field quantities describing compressible turbulence 
naturally generates a larger number of such correlations. These unknown correlations are 
all closed in terms of gradients of the mean fields and single point descriptors of the in- 
compressible field using Yoshizawa’s TSDIA formalism . 1 Unlike Yoshizawa’s theory 4 of 
compressible turbulence, which introduces the density variance as a new descriptor of 
strongly compressible turbulence, the present linearized theory only requires single point 
descriptors of incompressible turbulence. These are chosen following the usual two equation 
models as turbulence kinetic energy K and dissipation rate £. 

As Yoshizawa has demonstrated , 1 the TSDIA really generates an i nfin ite series for the 
unknown correlations both in powers of the mean gradients and in the order of differentia- 
tion of the mean fields. The present analysis is limited to the terms of lowest order in this 
expansion: those which are linear in the mean fields and which contain spatial derivatives 
of order two at most. A large number of potentially interesting nonlinear effects, analo- 
gous to the normal stress effects in simple incompressible shear flow, could be described 
by computing the higher order terms in the series. 
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In addition to the usual incompressible eddy viscosity, the calculation generates eddy 
diffusivities for entropy and pressure and an effective bulk viscosity acting on the mean 
flow. The compressible fluctuations also generate an effective turbulent mean pressure and 
corrections to the speed of sound. The “renormalized bulk viscosity, mean pressure, and 
sound speed are consequences of coupling among the response equations which correspond 
to the generation of sound waves in compressible turbulence. These effects would be 
suppressed by a “diagonal” approximation which ignores such couplings. 

These effective properties have the character of enhancements of molecular properties 
by turbulent fluctuations. But a prediction that appears to be unique to TSDIA is that 
terms containing gradients of turbulence quantities also appear in the mean flow equations. 
Because of the large number of such terms and their complexity, only a preliminary account 
of the effect of these terms on single point modeling is given. Explicit evaluation of these 
terms is reserved for future work. 

It should be stressed that the present model only attempts to close the correlations 
which arise from Reynolds averaging of the equations of motion. The effects of compress- 
ibility on the two-equation model 14 ' 15 are not addressed, although we indicate later how 
these effects could be analyzed whthin the present setting. 

An important question concerns the validity of linearized theories of compressible tur- 
bulence at asymptotically large times: 16 in the infinite time limit, the compressible field 
might build up to amplitude levels sufficient to invalidate a linearized theory. For exam- 
ple, Staroselski et al 8 argue that equipartition between compressible and incompressible 
fluctuations is a possible outcome of the long time evolution of compressible turbulence. 
But when the compressible field is generated entirely from the incompressible field, the 
weakness of the relevant energy transfer mechanism may permit a range of times, suffi- 
ciently large for many applications, over which the linearization remains valid. Numerical 
simulations of energy transfer in weakly compressible turbulence by Bataille et al 1 7 provide 
examples of this possibility. 

II. Formulation of the theory of weakly compressible turbulence 

The inviscid governing equations are taken in the form 

ds 

%- + u • Vs = 0 

dt 
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<9u 1 

— + u • Vu + -Vp = 0 

ot p 

dp _ 

— + u • Vp + 7 pV • u = 0 (1) 

with dependent variables entropy s, velocity u, and pressure p. Density can be evaluated 
from entropy and pressure through the thermodynamic relation 

s/Cv = log— - qlog— (2) 

Po po 

In what follows, 5 denotes entropy divided by Cv to simplify the notation. 

The analysis is based on a three-component decomposition of the velocity field into 
mean quantities, incompressible fluctuations, and small compressible perturbations. Ac- 
cordingly, first introduce the usual Reynolds decomposition of the fields s, u,p into mean 
and fluctuations and then further decompose the fluctuations into an incompressible field 
and compressible perturbations following Zank and Matthaeus. 2 Thus, the field quantities 
are written as 


s = S + s' 
u = U + u°° + u' 

p = P + p°° + p (3) 


where 5, U, P denote the mean entropy, velocity, and pressure and u 00 ^ 00 the incompress- 
ible velocity and pressure fluctuations. Note that the field U need not be solenoidal. Write 
the density as 

p = R + p' (4) 

The density fluctuation p' can be evaluated by linearizing Eq. (2) as 

p'=c-V 0 +p')-7~ 1 ^' (5) 


where the mean speed of sound is defined by the thermodynamic relation 



The compressible velocity perturbation u' need not be irrotational; in fact, the interaction 
of the compressible field with the base isotropic incompressible field u°° generates addi- 
tional solenoidal motions. The dynamics of the compressible perturbations s',u',p' will 
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be analyzed bv linearization about the mean fields and incompressible fluctuations, this 
defines the theory of weakly compressible turbulence. 

The mean equations are obtained by substituting the decomposition Eq. (3) into Eq. 
(1), dropping terms which are nonlinear in the compressible quantities, and averaging. 

Then 


— + U ■ VS = -V- < u°V > 

dt 

R(— + U • VU) + VP = -i?V- < u 00 !! 00 + u°°u' > -R< u 1 Vu°° > 
dt 

+ + p°V > -i < s'Vp™ > 

+ U • VP + 7 PV U = -V- < u°V° + u°y + qu'p 00 > 


dP 

dt 


+ (7 - 1) < u' ■ Vp°° > 


( 6 ) 


( 7 ) 


( 8 ) 


where in Eq. (7), p' has been eliminated in favor of p°°,p', and s' using Eq. (5). The 
derivation of Eq. (7) requires the expansion of the ratio 1 j{R + p ') in powers of p' / R. 
Following the usual practice in the theory of compressible turbulence, only the term linear 
in p' has been retained. Computing this series to higher order in p' would generate a 
series in powers of p°°/c 2 . Although linearization does not justify ignoring these terms, 
they generate a series in powers of the turbulent Mach number which can reasonably be 
truncated at lowest order. 

The incompressible field satisfies 

du °° + u ~ . Vu°° + -Vp°° = -[U • Vu°° + u°° • (VU) S ] + < u°° Vu°° > (9) 

dt R 


V ■ u°° = 0 


( 10 ) 


Eqs. (9) and (10) require that the incompressible pressure p°° satisfy the Poisson equation 

v • — Vp°° = -V • (u°° • Vu°°) - V • [U • Vu°° + u°° • (VU) S + < u°° • Vu°° >] 

R 

These equations have been investigated in Yoshizawa’s TSDIA analysis of incompressible 
turbulence; 1 the results of that analysis will be assumed in what follows. The traceless 
part of the mean velocity gradient appears in the incompressible momentum equation Eq. 
(9); it is defined by the decomposition 

w = (VP) S + (VP f 
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The component ( VU) C is added to the compressible momentum equation below. 

The derivation of the linearized equations for the compressible fluctuations is routine 
and leads to 


-x- + u°° • Vs' + C s = 0 
at 

+ u°° • Vu' + u' • Vu°° + -^Vp' - ^^Vp°° + - Ps'Vp 00 + C u 
dv' 

■— + u°° • Vp' + Rc 2 V • u' + u' • Vp°° + 7p°°V ■ u' + C p = -pT 


where the source term in Eq. (13) is the incompressible quantity 


+ dp* 

P T = -5- + u°° • Vp 00 

at 


and the terms coupling the fluctuating and mean fields in Eqs. (11 )-(13) are 


C 3 = U ■ Vs' + u°° • VS + u' • VS- < u°° • Vs' > 
C“=U-Vu' + u'-VU-^VP--^Vp'--L 


r p°°Vp c 


+ u°° • (VU) C - < u°° • Vu' + u'Vu°° - -^p'Vp 00 - jpp'Vp' > 
C p = U • Vp°° + U • Vp' + u°° • VP + 7p°°V • U + u' • VP 

+ 7 p'V • U- < u°° ■ Vp°° -|- u°° • Vp' + u' • Vp°° + 7 p°°V • u' > 


In Eq. (16), the density fluctuation must be expressed in terms of entropy and pressure 
by Eq. (5). 


III. DIA analysis of homogeneous isotropic weakly compressible turbulence 

Analytical theories of turbulence such as DIA 5 provide a systematic procedure for 
evaluating the correlations generated by Reynolds averaging. Although DIA can be for- 
mulated even for arbitrarily inhomogeneous and anisotropic turbulence, 18 the resulting 
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theory proves to be intractably complex. TSDIA provides a useful approximation. Heuris- 
tically, it computes perturbatively about an isotropic and locally homogeneous state of 
turbulence which is described to lowest order by DIA. The effects of coupling to mean 
fields are evaluated by a perturbation series in a scale ratio parameter. This perturbative 
treatment restricts TSDIA to weakly inhomogeneous and anisotropic turbulence, however, 
since universal single-point modeling of turbulence with arbitrarily strong inhomogene- 
ity and anisotropy is impossible in principle, this restriction does not unduly limit the 
applicability of TSDIA to turbulence modeling. 

Analytically, the local homogeneity of the lowest order TSDIA field means that mo- 
ments of this field at each point are evaluated using DIA results for homogeneous isotropic 
turbulence generalized for weak inhomogeneity by allowing all single point descriptors to 
vary slowly with space and time. This approximation, justified in the formal development 
of TSDIA, 1 occurs in many physical theories, for example in thermodynamics, when a 
thermodynamic system is taken to have a temporally or spatially variable temperature. 

In the present problem, Eqs. (11)-(13) describe a coupled system of passive fields. The 
application of DIA to this problem is not straightforward because the coupling between 
the compressible and incompressible fields cannot be analyzed exactly. Approximations 
are required which are both analytically tractable and preserve the important physical 
properties of the problem. Further discussion appears in Appendix I. 

The simplest suitable approximation is obtained by treating both the mean field cou- 
plings and the terms which depend on the incompressible pressure as perturbations. Thus, 
write Eqs. (11)-(13) as 

— + u°° • Vs = -C s (18) 

dt 

^ + u°° • Vu + u • Vu°° + — Vp = (p/pl)Vp°° - C“ (19) 

dt Po 

$ + u°° • Vp + pool V . u = -pt _ ( U - Vp°° + 7p°°V • u) — C p (20) 

at 

In Eqs. (18)-(20), the primes on compressible field quantities have been dropped. Since we 
are treating the homogeneous problem, R = po and c — cq are assumed constant in Eqs. 
(18)-(20). Following the perturbation scheme outlined in Appendix I, we must evaluate 



Note that the s equation Eq. (21) decouples from the others: s is a passive scalar 
advected by the incompressible velocity field u°°. The response matrix therefore has the 
structure 

' G aa 0 O' 

G = 0 G?f G up (24) 

0 Gi G pp 

The elements of the response matrix are defined as 5 

G ss =< Ss/Sf a > 

G =< Sm/SfJ > 

G u x p =< 6ui/6f p > 

G pu =< SpfSff > 

G pp =< 6p/Sf p > 

where the quantities /*,/*, / p are small perturbations added respectively to Eqs. (21)- 
(23). 

The spatial homogeneity of the fields justifies Fourier transformation in space. Since 
s is a passive scalar, its response function G sa satisfies the standard DIA equation 19 

dG 3s r f 

- 7 £-(k,f,r) + / ds ri aa (k,t,s)G aa (k,s,r) =6(t-r) (25) 
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( 26 ) 


where the damping factor p ss in Eq. (25) is evaluated as 

rj sa (k,t,r) = k m k n f dpdq G s9 {p,t,r)Q mn (qJ,r) 

J k=p+q 

In Eq. (26), Qmn is the incompressible two time correlation function 

< u^(k,t)u£°(k',r) > = Q mn (k, t, r)<5(k + k') 

The 6 function forcing in Eq. (25) summarizes the initial and causality conditions for G ss 
which are 

G SJ (M,r) = 0 for t < r 

G aa (M,f) = 1 

The spatial Fourier transforms of Eqs. (22) and (23) are 

^-(M) + — ;fcip(M)+ f dpdq T ! 1 mn (k, p, q)u^(p, t)u n {q, t) = 0 

ut p o J k=p+q 

?l(k,t) + ip 0 c 2 k i u i (k,t)+ l dpdqr 2 n (k,p,q)u~(pJ)p(q,t) = 0 (27) 

ut J k=p+q 

The couplings between the compressible and incompressible fields in Eq. (27) are defined 

by 

rj mn (k,p,q) = i(q m din +Pn6im ) 

T^(k,p,q) = ik n 

The DIA equations for the response functions are easily shown to be 

+ -ikG r * + v‘ P * Gf + r,** . G ;; = 6(t - r) 

3 po 3 

G“ p + — ik t G pp + T]* p * G pp + * G" p - 0 

Po 

G pu + ipoc 2 kjG U ji + p pu * G““ + V pp * Gf = 0 

G pp + ipoc 2 kjGY + v p j u * G up + n pp * G pp = S(t - r) (28) 

where * denotes time integration as in Eq. (25) and the damping functions t) are given by 


1 

II 

3* 

sr* 

/ dpdq r t 1 mn (k, p, q)G““ (q, t , r)rj r> (q, -p, k)Q mr (p, t, i 

J k=p+q 

pj p (k,t,r) = - 

f dpdq rj mn (k, p, q)G“ p (q, t, r)r 2 (q, — p, k)Q mr (p, f, r) 

J k=p+q 

p pu (k ,t,r) = - 

f dpdq r 2 m (k, p, q)G pu (q, t, r)rj, ri (q, -p, k)Q rm (p, t, r) 

J 

J k=pH-q 

Tj pp (k, t.r) - — 

f dpdq r^(k, p, q)G pp (q, t,r)r^(q, — p, k)Q mn (p, f , r) 

Jk=p+q 
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Now introduce the normalized Helmholtz decomposition of the field u*: 

ti,(k, <) = u>,(k, t) + ikik~ 1 4>(k, t) 

The isotropy of the lowest order fluctuating fields implies that the damping and response 
functions must have the form 

7 ? ! U /(k. t, r ) = r i ww (k, t, r)Pij( k) + t, r)P* } { k) 

r]f p (k,t,r) — rj <l>p (k,t,r)ikik ~ 1 

T] pu (k ,t,r) = -r) p4, {k,t : r)ikjk~ l (30) 

and 

k, *, 5 ) = G™ a)P 0 - (k) + G++(k, *, s)P* } ( k) 

G“ p (k,M) = G* p (k,t,s)ikik~ l 

G pu ( k, f, a) = -G p *(£, f, a)**,-*- 1 (31) 

In Eqs. (30) and (31), P tJ and P* are the transverse and longitudinal projection operators 

-P,y(k) — Sj \j kikjk 
P*j(k) = kikjk - 2 

Substitute Eqs. (30) and (31) into the governing equations Eq. (28) and separate 
transverse and longitudinal components. The result is that G ww decouples from the rest 
of the system and satisfies, in the abridged notation of Eq. (28), 

G ww + vf ww * G ww = 6{t - r) (32) 

while the remaining response functions satisfy the system 


J7 1 Ol d [ 0 G* p 

{ [o l\dt + [-poc 2 k 0 Jj[G^ G pp 

rj+P 1 G** G 4,p ]_ fl O' 

+ [rp+ t]Pp\ * [G p4> G pp J _< ^ r) [o 1 


(33) 


10 



The solution of Eq. (33) for the response equations is discussed later. For now, the response 
matrix with the simplified structure of Eq. (24) is used to derive formal expressions for 
the transport coefficients. 


IV. TSDIA analysis of inhomogeneous weakly compressible turbulence 


TSDIA introduces, in addition to the usual Reynolds averaging, a two-scale decom- 
position 


d_ _ d_ _d_ 

dxi dx[ dXi 

d_ _ d_ d_ 

dt~ dt' + dT 


(34) 


where x' and ? describe small turbulent scales of motion and X and T describe large scale 
flow features including inhomogeneity and anisotropy. It is natural to assume that the mean 
quantities depend only on the slow variables X, T. Substitution of the scale decomposition 
Eq. (34) in the mean equations Eqs. (6)-(8) therefore leaves these equations unaltered. 
Corresponding to Eq. (34), write the compressible field quantities as 


s = s(x',X,f',T) 
U{ = Ui {x',X,t',T) 
p = p(x',X,t',T) 


and apply the scale decompositions of Eq. (34) to the equations for fluctuations Eqs. (18)- 
(20). With these substitutions, dropping primes in the small scale variables, the fluctuation 
equations take the form 

— + V ■X x s + u 00 -V x s = -r (35) 

at 

— + U • V*u + u°° V.u + u • V lU ~ + lv,p = -r- (36) 

di n 

^ + U- V lP + u°°- V xP + Rc 2 X x -u = -X p (37) 


where 



F u = {-QT + U ■ Vxu) + u°° ■ V,vu + u V.vu°° + jjVxp 

- w v - p0C - W Vxp °° + u • v * u - w VxP - ~k Vxp - F v * p 

+ u *-(V x Uf (39) 

T* = (jfi + U • Vxp) + u°° • V x p + Rc 2 V x ■ u 

+ u • V.y p°° + u • X x p°° + 7P°°^A' ■ u + 7 p°°V x • u 

+ + U ' Vxp °°} + U ‘ V ^°° + u °° • VvP + 1P°°V X ■ U 

dn°° 

+ 7 pV A - • U + u - V X P + u°° • Vxp°° + (-^- + u°° • V^p 00 ) (40) 

In Eqs. (35)-(40), the variable mean density R and sound speed c appear because we now 
consider weakly inhomogeneous fields. The averages in Eqs. (15)-(17), which maintain 
the property that the fluctuation equations have zero mean, have been dropped in Eqs. 
(38)-(40) because they do not contribute to the formation of moments at the level of 
approximation to be introduced later. 

Since the lowest order field is locally homogeneous, a spatial Fourier transform is again 
appropriate. The spatial part of the decomposition Eq. (34) becomes 

d d 

dxi * dXi 

With these assumptions, the fluctuating fields are expressed in terms of the corresponding 
forces and response functions as 


s(k,X,#,T) = - / dr G 3S (k,X,t,r,T)^ s (k,X,r,T) 
Jo 

(41) 

Ui (k,X,LT) = - f dr [G“ u (k,X,f,r,T)^“(k,X,r,T) 


+ G“ p (k,X,t,r,T)^(k,X,r,T)] 

(42) 

p(k, X, t, T) = - [ dr [G pp ( k, X, t , r, T)T P ( k, X, r, T ) 
Jo 


+ Gf(k,X,t,r,T)J7(k,X,r,T)] 

(43) 


The possibility of writing explicit expressions like Eqs. (41)-(43) is an important advantage 
of the DIA formalism. But it must be emphasized that these expressions can only be used 
to form moments containing s,u*, and p: DIA provides statistical formulas for the field 
quantities, not pointwise solutions of the equations of motion. 
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Analytically, local homogeneity means that the various response and correlation func- 
tions G(k, X,t, r, T) and Q(k,X,t,r,T) are evaluated as the corresponding homogeneous 
quantities G(M,r) and <?(M,r) but with all single point moments which characterize 
the fields treated as functions of A', T. For example, consider the typical forms 1 for the 
response and correlation functions of isotropic incompressible turbulence, 

G{k,t,r) = exp [~T}(k)(t - r)\H(t - r) 

Q{k, t, r) = ^-C K £ 2 / 3 k~ n/ 3 exp[-r](k) \ t - r |] 

47T 

v (k) = C D e^k^ 

k>k 0 ( 44 ) 

where if is the unit step function. These expressions are generalized to weakly inhomoge- 
neous incompressible turbulence by letting £ and ko depend on X and T. 

With these results, it is straightforward to evaluate the moments needed to close the 
mean flow equations. In evaluating the moments, the following truncations of the series 
which the TSDIA generates are adopted: 

(a) only modeled terms which are at most linear in the mean gradients are retained, 

(b) only modeled terms with at most second order derivatives with respect to A', T are 
retained, 

(c) modeled terms derived by repeated convolution with the compressible response func- 
tions are dropped. 

Whereas (a) and (b) follow from the formal development of TSDIA as a perturbation 
expansion* (c) is an additional approximation. 

In the following, dependence on A', T will be understood but not written explicitly. 

A. Entropy velocity correlation 

The entropy velocity correlation appears in the mean entropy equation Eq. (6). Write 
this correlation as the integral 

< > = J dk < ur(-MMM) > (45) 
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Substitute the expression for s in terms of its response function Eq. (41) in Eq. (45) to 
obtain 

-<uf°s> = J dk dr G S3 (k,t,r)x 


r) o 

{<u~(-M)^(k,r)> 

(46a) 

ds 

+ <u°°(-k,t)C/ p — (k,r)> 

(466) 

+ < u“(-k,t)u“(k,r) > 

(46c) 

+ <«r(-M)«j(k,r)>|l} 

(4 6d) 


The formation of nonzero correlations between u°° and compressible fields in lines 
(46a), (466), (4 6d) requires an additional time integration over the response functions; the 
approximation (c) above suppresses these terms. It can be noted that some of these terms 
will also generate products of gradients and gradients of higher order which are also ex- 
cluded by (a) and (b). 

Term (46c) contributes a gradient transport expression 

- <u?(-kj)s(kj)> = [ dk / dr G a ’(kXr)Qij(k,t,r)4£r (47) 

J J o uX j 

To lowest order of TSDIA, the integral in Eq. (47) is evaluated by substituting isotropic 
quasi-homogeneous forms for G 3S and Q t] : 

Qii( k,t,r)=Q(k,t,r)P i] (k) 

G 33 (k,t,r) = G 33 (k 7 t,r ) (48) 

To construct a single point model, it is assumed that the fast time scale dependence 
is stationary. Therefore, the response and correlation functions depend only the time 
difference t = t — r: 

Q{k,t,r) = Q(k,t - r) = Q(k,r) 

G ss ( k. t, r) = G 33 (k , t-r) = G 33 (k , r) (49) 
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Furthermore, since the turbulent time scales are much shorter than those of the mean 
motion, the long time limit t -* oo is taken. Substituting Eqs. (48) and (49) in Eq. (47) 
and taking the long time limit results in the gradient transport model 


- < urs > - v 


, dS 
dX t 


(50) 


Avhere the transport coefficient is defined by 


V 


ss 



47 rk 2 dk 



dr G sa (fc,r)<2(fc,T) 


(51) 


The factor of 2/3 arises from angular integration of the transverse projection operator 
Pij. This type of expression for a transport coefficient is familiar from Yoshizawa’s TSDIA 
calculations. 1,4 


B. Velocity pressure correlations 

Four such correlations appear in Eq. (8) for the mean pressure: < u 00 p 00 >, 
< u°°p >, < up 00 >, and < u Vp°° >. 


1. The correlation < u 00 ^ 00 > 


The incompressible correlation < u 00 ]? 00 > has been analyzed by Shimomura 20 who 
finds that it is a diffusion of turbulence quantities. In terms of single point quantities, the 
integrals become 

(52) 


< P °°u? 


A' 2 dl< , c Kde_ ) 


ax. 


Terms of this type will also appear in the closure of correlations containing compressible 
fields and will be discussed below. 


2. The correlation < u°°p > 

Write the second velocity pressure correlation as the integral 

- < u°°p > = - J dk < u°°(~M)P(M) > 

= J dk j* dr {G pp (k.t,r) < u?{-kJ)T p (k,r) > 

+ G pu (k.t,r) < u?(-k,t)fj(k,r) >} 


(53) 
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Closure of this moment, requires that the forces from Eqs. (38)-(40) be substituted in Eq. 
(53). The approximation in which only one time convolution with a compressible response 
function is allowed shows, as for entropy correlations, that we need only consider terms 
containing correlations of the incompressible field with itself; the result of evaluating only 
such terms appears in Appendix II, Sect. A. 

Further analysis shows that some of these terms vanish. A list of the nonvanishing 
terms also appears in Appendix II, Sect. A. Among these terms, some have the form 
familiar from eddy viscosity modeling: they are (1) proportional to gradients of mean fields, 
and (2) nonvanishing in regions of constant turbulence properties. Property (2) gives these 
terms the character of enhancements of molecular properties due to turbulent fluctuations; 
the molecular properties themselves are the result instead of thermal fluctuations. The 
transport coefficient v ss of Eq. (51) is an example. Such terms will be denoted by the 
subscript M; only these terms will be evaluated in detail in this paper. 

But TSDIA also generates a second set of terms which depend on gradients of tur- 
bulence quantities. Shimomura's closure of the pressure velocity correlation Eq. (52) is 
of this type. Unlike the terms just described, these terms vanish when the turbulence 
properties are constant. Since turbulence properties are rarely constant in practice, it is 
an important and nontrivial prediction of TSDIA that such corrections to the transport 
model exist. Because of the large number of such terms and the lengthy analysis required 
to evaluate many of them, we defer a complete analysis of these terms to later work. In 
this paper, such terms will be labelled with the subscript T. A preliminary description of 
them appears in Appendix III. 

With these conventions, 


< U p > = < U p > M + < U °°p > T 


It is shown in Appendix II, Sect. A that 


- < u°°p >m= v pp 


dP 

dX t 


(54) 


where, repeating the argument which led to Eq. (51), the transport coefficient v p x p is given 
by 

2 /*oo /»oo 

V \ V = 3 y o 47 rk2dk l dr G pp (k,r)Q(k,r) (55) 
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3. The correlation < up 00 > 


Write the third velocity pressure correlation as the integral 

- < mp 00 > = - j dk < p°°(-k,t)ui(k,t) > 

= j dk f dr {GHM,r) <p~(— k,t)**(k,r) > 

+ G“/(k,t,r) <p°°(-k,^;(k,r) >} (56) 

The terms generated by substituting the forces of Eqs. (38)-(40) in Eq. (56) are listed in 
Appendix II, Sect. B. The result has the form 

< up 00 > = < up°° >Af T < up°° >T 


where the mean flow dependent term is 


— < u lP °° >m — 



dP 

dXi 


(57) 


and the transport coefficient v pp is given by 

^ = J_/ 47 rk 2 dk dT{^G ww (k,T) + -G^(k,T)}Q p (k,r) (58) 

R 2 C 2 Jo Jo 3 J 


In Eq. (58), we have assumed, corresponding to Eq. (49) for the two-time velocity corre- 
lation function, an isotropic stationary form for the two-time pressure correlation 


< p°°(k, t)p°°( k, r) > = Q p (k, t)S ( k + k') 

Whereas the single time pressure correlation function is described by Batchelor’s analysis, 21 
calculation of the two-time pressure correlation remains an open problem. 

Eq. (57) has an anti-diffusive character, but the complete transport coefficient is 
v pp = - vl p where v\ p is evaluated in Eq. (55); it will be shown later that the negative 

contribution is smaller in magnitude than the positive contribution by a factor of the order 
of the turbulent Mach number squared. 


17 



4 . The correlation < u ■ Vp°° > 


The fourth pressure velocity correlation is written as an integral 


< u • Yp°° > 


■ / dk a 


dr 


{ < GHM.r).P’(k,r)(-<* j + j|-)p“(-k,i) > 

+ < G?;(k,(,r)^7(k,r)(-.'ty + g|-)p“(-k,<) >} 

It is shown in Appendix II, Sect. C that the only mean flow dependent term is 

dU m 


< u v x p°° >m = n 


dX , 


where the effective pressure II is given by 

r oo roc 

n = 7 / 47rPdA: / dr kG^ik, T)Q?(k, r) 
Jo Jo 


( 59 ) 


( 60 ) 


C. Velocity-velocity correlations 

The velocity-velocity correlations arise in the mean momentum equation Eq. (7). 
There are three such correlations: < u^u 00 >, < u°°u >, and < u • Vu°° > 

1. The correlation < u 00 !! 00 > 


The incompressible correlation < u^u 00 > has been extensively investigated by 
Yoshizawa. 1 At the level of approximation evaluated here, the result of the TSDIA analysis 
is the usual linear eddy viscosity representation 


- < t*r«r > = -\k6h + v(s l} ~ 


( 61 ) 


where Sij is the strain rate 


C . _ dU > QU J 


11 dX } ' dXj 
and the incompressible eddy viscosity is the integral 1 

4 r°° r°° 

v = — / Airk 2 dk I dr Q(k,r)G(k,T ) 
15 Jo Jo 

where G is the incompressible velocity response function. 


( 62 ) 
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2. The correlation < u°°u > 

The second correlation is written as the integral 


- < ufuj > 



< uf°(-M)u,-(M) > 
dr {G“"(k,*,r) < <°(-k,*)JF“(k, 


r) 


+ G“ p (M,r) < u?(-k,t)F p (k,r) >} 


> 


Reference to Appendix II Sect. D shows that there is only one mean flow dependent term, 
which can be described as a renormalized bulk viscosity 


— < u t Uj > m =Vi 


UU ^ 


dX, 


(63) 


where the transport coefficient is given by 


v 


UU 

1 


1 f°° 

- / 4t rk 2 dk 

3 J 0 



dT^G WW (k,T)Q(k,T) 


(64) 


3. The correlation < u • Vu°° > 

This correlation is written as the integral 

f d 

- < (u • Vu°°), > = - I dk < «j(M)(-*'*j + — )wr(- k ^) > 

= j dk J dr {G**(k,f,r) < JZ(k,r)(-,kj + ^KVM) > 

+ Gj p (k, t,r) < X p (k ,r)(-ikj + — )u°°(-k,t) >} 

It is shown in Appendix II. Sect. E that the gradient transport term is 

_<(u.Vu (65) 

where the coefficient R is defined by 

j* oo roo 

4= / An k 2 dk / dr fcG* p (fc,r)Q(fc,r) (66) 

R Jo Jo 

It will be shown later that R is related to an effective or renormalized speed of sound waves 
propagating in the mean flow. This effect has been predicted b\ Chandrasekhar and b\ 
Staroselski et al. 8 
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D. Pressure density correlation 

There are three terms: < p^p 00 >, < pp°° > and < sVp 00 >. 

1. The correlation < p 00 ^ 00 > 

The incompressible pressure correlation < p^p 00 > is the subject of Batchelor’s classic 
analysis 21 which gives 

< > = CpK 2 (67) 

where the constant Cp can be evaluated theoretically or from measured data. This result 
also holds in the lowest order of TSDIA. 


2. The correlation < pp°° > 

For this correlation, 

— < p°°p > = — J dk < p°°(— k,t)p(k, t) > 

= J dk I dr{G pp (k,t,r) <p c °(-k,t)r'(k,r) > 
+ G pu (k.t,r)<p°°(-k,t)r’(k,r) >} 

In Appendix II Sect. F, it is shown that the mean flow dependent term is 

- < pp°° >m = R 2 c\r^f- 

where the effective bulk viscosity is defined by 

JJ 2 C 2 I/J“ = 1 J i-k-tll' J dr G pp (k,T)Q p (k,T) 


( 68 ) 


(69) 


Note finally that since there is no pressure source in F s , the correlation < sp°° > 
vanishes to this order. 

Introducing only the mean flow dependent terms, Eqs. (50), (54), (57), (59), (61), (63), 
(65), and (68) computed by TSDIA into Eqs. (6)-(8), we obtain the gradient transport 
model for weakly compressible turbulence, 

A 

— + U . V5 = V • (I/-V5) (70) 

^ + u-yu + (i- i)vp = v • {-|at + i/[vu + (vu) T ] 5 } + V(V uu v • uni) 

^ + U • VP + ( 7 - 7)PV • U = V • (u pp VP) (72) 
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The gradients of turbulence quantities generated by TSDIA can of course be added as well. 
A preliminary qualitative description of these terms and a discussion of the form of the 
resulting model appears in Appendix III. 

The transport coefficients in Eqs. (70)-(72) are found from Eqs. (51), (55), (58), (60), 
(62), (64), (66), and (69) above and are summarized here for convenient reference as 


a r°° r°° 

U — — / 47T k 2 dk / dr G(k,T)Q(k,T) (73) 

15 Jo Jo 

cy r oc roc 

u a3 = - 4:irk 2 dk dr G S9 (k,T)Q(k,T) (74) 

3 Jo Jo 

9 r oo r oo 

v pp = - 4tt k 2 dk drG pp (k,T )Q(k, t ) 

3 Jo Jo 

i r°° o 1 

_i/ 4 7rk 2 dk dr{-G^(k,r) + -G^(k,r)}Q p (k,r) (75) 

C 2 Jo Jo 3 J 

7 = (7 - l)Xl/P < 76 » 

n = 7 / 4 tt k 2 dk / dr kG 4 ’ p {k,T)Q p {k,T) (77) 

Jo Jo 

U uu = - 47 rk 2 dk / dr^G ww (k, r )<5(^, r ) ( 7S ) 

3 Jo Jo 3 

_ 7 / 47f k 2 dk / dr G pp (k,r)Q p (k,r) 

R 2 C 2 Jo Jo 

1 TOO 

A= Uk 2 dk dr kG* p {k,T)Q{k,r) (79) 

R Jo Jo 


The eddy viscosity ^ and eddy diffusivity are transport terms of the type familiar 
from studies of incompressible turbulence. The term u uu is an effective bulk viscosity- 
generated by the compressible fluctuations. The modified pressure II has been introduced 
into the transport model Eq. (71) through the modified specific heat ratio 7 defined by 
Eq. (76). The turbulent fluctuations therefore modify not only the viscosities seen by the 
mean flow, but also the effective 7. Of course, the thermodynamic specific heat ratio is 
unchanged. The TSDIA expressions for II demonstrate that this modification of 7 is a 
long time effect: at short times, II ~ t 2 , hence 7 ~ 7 + 0(t 2 ). 

To clarifv the role of the R term, note that if the mean field is treated as an ensemble 
mean, then it can contain non-random sound waves. The propagation speed of these waves 
can be found by linearizing the mean flow equations about a constant state, setting 

P = Po + P' 
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U = Uo + U' 
R = Rq -f- R 1 


where Pq and R 0 satisfy the thermodynamic relation c 0 = (jPo/Ro) 1 ^ 2 where c 0 is the 
sound speed at constant pressure P 0 and density R 0 . The standard derivation shows that 
sound waves propagate in the mean field with speed 

c 2 = [7P0 - (7 - l>n][i- - i] = c 2 [l - 7(7 - 1) JL][1 - *L] (80) 

By evaluating the terms II and R explicitly later, it will be shown that Eq. (80) corresponds 
to an increase in the effective speed of sound in the mean flow. An argument for the en- 
hancement of the speed of sound by compressible fluctuations is given by Chandrasekhar; 22 
Staroselski et al 8 discuss this effect in the case of strong compressibility. 

V. Evaluation of the transport coefficients 

Explicit expressions for the response functions are needed in order to evaluate Eqs. 
(73)-(79) for the transport coefficients. In TSDIA calculations 1 for incompressible turbu- 
lence, the required expressions are given the Kolmogorov inertial range forms Eq. (44). In 
the present case, it is necessary to solve Eqs. (28) for the compressible response functions. 
The solution of these equations is difficult because the damping factors rj of Eq. (29) 
depend on the response functions themselves. Although these equations are already the 
result of several approximations, the desirability of analytical expressions for the transport 
coefficients justifies further approximations. 

First, it will be convenient to replace the damping by Markovian damping by setting 

r]{k,t,r) = r}(k)8(t - r) 


For time stationarity. appropriate to the lowest order TSDIA field, 

r OC 

rj(k) = dr rj(k,T) 

Jo 

where r = t — r denotes time difference. Equivalently, in frequency space 

T](k) = T](k,0) = rj(k) 


( 81 ) 
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This Markovianization is familiar in statistical physics 23 and is implicit in all turbulence 
modeling. It essentially states that turbulent processes occur much faster than mean flow 
processes. When this assumption is not valid, universal single time turbulence modeling is 
impossible in principle. Further discussion can be found in the recent work of Yoshizawa 24 
and Woodruff. 25 We note that the fluctuation-dissipation relation which connects the two- 
time correlation and the response function in Eq. (44) is justified for Markovianized 
theories. 25 

The time stationary form of the response equations Eq. (33) which includes the initial 
conditions is 



—iuj k/po 


e- 

•S- 

^3 

\ 

’ q4>4> 

G™' 


'l o' 


—poc 2 k —iu) \ 

+ 

^3 

0- 

1 

/ 

QV<t> 

Gpp 


0 1 


The Markovianized damping factors Eq. (29) with Eq. (81) can be rewritten as 
V ww (k) = k) / dr f dpdqrJ mn (k,p,q)G u;u ’(q^)x 

~ Jo J k=p+q 

-Pnp(Q)r^ r j(q 7 ”P? k)Q mr (p, t ) 

-l-Piji k) / dr f dpdqr i 1 rnn (k,p,q)G^(q,r)x 

2 Jo J k=p-fq 

P^(q)rJ rj (q, -p,k)Q mr (p, r) 

v **(k) = -P*(k) [ dr f dpdqT] mn (k,p,q)G ww (q,T)x 

Jo J k=p+q 

P np {q)T), rj (q, -p. k)Qmr (P, r) 

r oo r 

-P*(k) / dr dpdqrLn(k,P,q)<? <W, (q,T) x 

J Jo J k=p+q 

P* p (q)rJ rj (q, -p, k)Q m r(p, r) 

rj 4 >p {k) = ik r k ~ 1 f dr / dpdq rj mn (k, p, q)G“ p (q,r)r 2 (q, -p, k)Qmr(p, t) 

Jo J k=p+q 

r 7 P<fr (k,r) = — ikjk~ l f dr j dpdq T^(k, p, q)G pa (q, , r)rJ irj (q, -p, k)Q rrn (p, r) 

Jo j k=p+q 

r/ pp (k, r) = - [ dpdq T^(k, p, q)G pp (q, r)r 2 (q, -p, k)Q mn (p, r ) (83) 

J k=p+q 

The Markovianized system Eq. (82) for the response functions is solved approximately 
by perturbing about purely incompressible damping. Note first that Eq. (32) for G uu 
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reduces to the equation for the response function of a passive vector field if time integrals 
containing G 4>4> are ignored in comparison to time integrals containing G ww in the integral 
for rj ww in Eq. (83). The self-consistency of this approximation will be demonstrated later: 
it will be shown that a consequence of this approximation is that time integrals in G 4>4> are 
of order c~ 2 times time integrals in G ww . So approximate G ww by the response function 
of a passive vector field 


G ww (k , f, s) = exp [-T] ww (k)(t - s)}H(t - s) 


(84) 


where the damping factor r/ ww 


is related to the incompressible damping function rj through 


T] ww (k) = a^k) 

r)(k) = C D £ 1/3 k 2 ^ (85) 


In Eq. (85), a w is the inverse Prandtl number for convection of w. This quantity can 
be evaluated theoretically following the calculation of the Prandtl number for passive 
fields. 26,27 Eqs. (83) show that the integrals for the damping factors r ) <t>p , r ) p<t> , rj pp all 
contain the corresponding response functions G^ p , G p4> , G pp . Only the integral for tj 4 ’ 4 ’ has 
a contribution which depends on G ww . Assume further that time integrals which contain 
the response functions G <pp . G pd> , G pp are small compared to time integrals which contain 
G ww . Then set as a first approximation 

»?^( k ) = a^rjik) 

71* P (k) = 7 f+(k) = T] PP (k) = 0 

where a $ is another inverse Prandtl number which can be computed by evaluating the 
integrals in Eq. (83). Like Eq. (84), it will be shown subsequently that this approximation 
is self consistent: by computing response functions and corrected damping factors, we will 
find that if p , r} p<t> ,r] pp are of the order c~ l times 7]^. With this initial approximation, the 
solution of Eq. (82) is 


G* p {k) = 
G pp (k) = 
G*+(k) = 
G p4> {k) = 


Po 1 ^ 

— ui 2 + c 2 k 2 — iu>a ( f > r}(k ) 

—iu) + o^r]{k) 

—uj 2 + c 2 k 2 — iuia^r]{k) 

—ito 

—u 2 + c 2 k 2 — iida^ri(k) 

PoC 2 k 

— id 2 + c 2 k 2 — iida ( f,ri(k) 


(86) 
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The corresponding time domain expressions for the response functions are 
G^ik.r) = exp(— o;^r;(A:)r/2)sin(5r)ff(r) 

p 0 S 

G p4, (k,T ) = po —^-Gxp(—a ( j > r)(k)T /2)sin( Sr)ff(r) 

G pp (k, r) = exp(— a^ikjr /2 )cos(St )H{t ) 

+ 2^exp(-a^(t)r/2)<in(ST)J(r) 

G <t><t ’{k , r) = exp (— a ^r] (k )t / 2)cos( St )H(t ) 

_ a^) eX p(_a^(A:)r/2)sin(5r)F(T) (87) 

where 

s = { C H* - (88) 

This approximation can be improved by substituting the response functions of Eq. (8 1 ) into 
the formulas Eq. (83) for the damping factors. The corrected values of r^ p , t] p4> , r) pp contain 
terms which oscillate rapidly with frequency ck\ time integrals of these quantities will 
consequently be of order c" 1 . Consequently, the corrected i/*', V pp are of the order c" 1 
times r}^ . By solving Eq. (82) for the response functions with this second approximation 
for the 77's, a second approximation to the response matrix is generated. This iteration 
generates a sequence of approximate solutions to Eq. (82) containing increasingly high 
order powers of the time scale ratio r)(k)/ck. The approximation Eqs. (86)-(87) is the 
approximation of lowest order. This approximation has also been stated by Bertoglio et 

al. 6 

With these approximations, it is possible to evaluate both the transport coefficients 
of Eqs. (73)-(79) and the correlations of compressible field quantities which enter the 
two-equation model for weakly compressible turbulence. These correlations, such as the 
pressure-dilatation and dilatational dissipation, are believed to be crucial to explaining 
compressibility effects on mixing layer growth. To evaluate them, it is necessary to solve 
Eqs. (18)-(20), with all mean field contributions dropped, to express the compressible 
fields as integrals of the response functions of Eq. (86) and the incompressible source term 
pt of Eq. (14). The result will express these compressible correlations in terms of the 

two-time correlation 

Qt.(k,t,r)$(k + k') = < pt(k,<)pt(k',r) > 
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Although this program is completely routine, Q t is a sixth order correlation of velocity 
components, therefore explicit evaluation of the required correlations will be quite lengthy. 
This part of the model development is left to subsequent work. 

To evaluate the transport coefficients, we apply the inertial range formula Eq. (44) 
with the single time correlation function given by the cutoff Kolmogorov spectrum 

E(k) = Q(k)HTk 2 =\CK eVH ~ SI3 , (89) 

( 0 if k < k 0 or k > kj 

In Eq. (89), kd is proportional to the Kolmogorov scale; the convergence of all integrals 
in this theory as k d — +• oo implies that we can set kd = oo and ignore its finite value. The 
large scale cutoff k 0 is defined in terms of the kinetic energy of turbulence K through 

K = ^C K e 2/3 k~ 2/3 

We take also 

G ss (k,r ) = exp[-arf(k)T]H(T) ( 90 ) 

where a is the inverse Prandtl number for a passive scalar. 27 Pending a more complete 
investigation, we postulate for the two time pressure correlation 

Q P ( k - T ) = -^^E p {k)exp[-apr]{k) | r |] (91) 

where the single time pressure spectrum is given by Batchelor’s well-known relation 

E p (k) = C B pie 4 '*k-V* (92) 

The spectrum constant C# and Prandtl number 1/ap are evaluated theoretically or from 
data. In evaluating integrals dependent on the term 5 of Eq. (88), we will make the 
leading order approximation 

S ~ ck 

which is consistent with the assumption of small turbulent Mach number. Given these 
approximations, it is easy to evaluate the transport coefficients in terms of single point 
quantities. 
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The coefficient v ss of Eq. (74) is evaluated as 


cy p OC pOO 

= - / dk dr G ss (k,T)E(k)exp(-riT) 

3 Jo ./o 

— f 

+ 1 Jo 




2 

I dk ... 

3 a + 1 J 0 r](k) 

2 1 Ca~ 1/ 3 ,-4/3 

3 a + 1 0 

9 Q + 


— (Ca-Cd) 1 — 

1 £ 


(93) 


The transport coefficient i/ ss describes the convection of a passive scalar, and is therefore 
of order M t ° where the turbulent Mach number is M t - I< 1/2 jc. 

The coefficient u pv is evaluated from Eq. (73) as 


v 


■pp 


cy pOO pOO 

= - dk drG pp (k,T)E(k)exp(—rfT ) 

3 Jo Jo 

1 roc poo cy 1 

__L/ dk dT{^G ww (k,T) + -G <t> *(k,T)}E p (k)exp(-apr } T) 

C z Jo Jo 3 3 

f°° 2 77 1 2 1 E p (k) 

= J 0 ^3 ( ^ + 1) ?F £(J:) "?3ap + «„ »7(fc) 

= [i( Q . + l)C K C D -l-i-§|]^- 


= |-c? 3 [(«* + 1)C K C D - 


1 


— — 

&p + oi w Id £ 


-2 


(94) 


This evaluation is correct to lowest order in M<; including the contribution from would, 
in view of Eq. (86), lead to corrections of order M 4 . Diffusion of pressure is therefore an 
O(Mf) phenomenon. 

In evaluating II from Eq. (77), we will again compute to leading order in M t . Thus, 

pOO poo 

n= 7 / 4irk 3 dk / dr G^(fc,r)Q p (fc,r) 

Jo Jo 

= - 7 / 

Jo 


oo rt 

dA: J^e 4 /^" 7 / 3 + 0(A/ 4 ) 


= -7‘ 


2poc 2 

< (P°°) 2 > 
2poc 2 


(95) 


where Eq. (67) can be substituted for the total pressure variance so that 


n ~ -KM 2 


(96) 
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For the bulk viscosity v uu , Eq. (77) implies 


v 


uu 



oo 

4nk 2 dk 


4nk 2 dk 


l 


/•oo 9 

J dT-G"“(k,T)Q(k, T) 

o 

cfr G pp (k,T)QP(k,T ) 


( 97 ) 


Note that the negative contribution to Eq. (97) is of order c~ 2 compared to the first, 
therefore it suffices to evaluate the first term only. The result is 


v 


UU 



dk 2 1 E(k) 

9 a w + 1 r](k) 


1 1 Us , -A/3 

6a w + lC D 0 


2 1 

2 ( ot w + 1 


(CkCd) 1 



(98) 


Note the somewhat surprising conclusion that the bulk viscosity has an 0(M < ° ) contribution 
from the fluctuating w field. The reason is that the decorrelation time of the w field appears 
in the expressions for the transport coefficients, not its amplitude, and that in the present 
approximation for the response functions, w decorrelates like a passive vector. Although 
the bulk viscosity does not vanish in the limit M t — > 0, incompressibility is maintained 
through the condition that V • U — > 0. 

Evaluating R to leading order in M t from Eq. (79), 


1 

a 



4:Trk 2 dk 



dr kG <l,p (k,T)Q(k,T ) 


1 f°° 

-- dk c~ 2 E(k) 

R Jo 

- 4 m < 2 

R * 


(99) 


Combining the results of Eqs. (96) and (99) in the formula Eq. (80) for the effective speed 
of mean flow sound waves, 

C 2 = c 2 [ 1 + M 2 ] 

since Tl/c 2 ~ M? according to Eq. (96). This result agrees up to the constant multiply- 
ing M 2 with the prediction of Chandrasekhar 22 for the speed of propagation of density 
correlations. 


28 



To summarize, by substituting the leading order terms in M t from Eqs. (93)-(99) for 
the transport coefficients in Eqs. (d3)-(t2), we obtain the mean flow equations 

<99 K 2 

+ U • V5 = V • (C s — VS) 
dt £ 

— + U ■ VU + i[l + Mt]VP = V ■ {-?A'I + C,-[VU + (VUff) 

dt R 6 s 

+ V(C U ^-V-U) ( 10 °) 

— + U • VP + 7[1 + C 7 M t 4 ]PV • U = V • (C P M 2 VP) 

dt 6 

where the model constants C s , C 7 , C p , CV, C u are given by explicit integral expressions 
containing the velocity inertial range constants C K and C D and the incompressible pressure 
inertial range constants Cp and ap. 

VI. Improved approximate system of equations 

Eqs. (21)-(23) were derived by treating all incompressible pressure terms as perturba- 
tions. This approximation lead to the decoupling of the entropy equation and consequently 
to the simple partitioned structure Eq. (24) of the response matrix. If instead all of the 
pressure terms in Eqs. (18)-(20) were treated non-perturbatively, the appearance of en- 
tropy in the momentum equation would cause the response matrix computed by DIA to 
be a completely filled three by three block matrix. Even an approximate solution of the 
corresponding response equations would be more difficult. 

However, the decoupling of entropy and the convenient matrix structure of Eq. (24) 
can be retained by treating only the term coupling s and p°° as a perturbation. This 
approximation is particularly reasonable since it is found in numerical simulations 28 that 
homogeneous compressible turbulence is nearly isentropic. With this approximation, the 
response equations still have the form of Eq. (28), but the damping factors become 

r;“"(k, t, r) = - f dpdq r j 1 mn (k,p,q)G“p(q,t,r)rp f . J (q,-p,k)Q mr (p,<,r) 

J k = p+q 

+ -~r [ dpdqp l k ] Q p {p,i,r)G pp (q,t,r) 
c ~ R A=p+q 

r/“ p (k, t, r) = - f dpdq r) mn (k, p, q)G“ p (q,t, r)T 2 p (q, -p, k)Qm P (p, t, r) 

J k=p+q 
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+ ~T^2 f d P d qp,PjG pu (q,t,r)Q p (p,t,r) 

c rL Jk=p+c\ 

T! pu (k,t,r) = - j dpdq T 2 p (k,p,q)G pu (q,t,r)T 1 nrj (q > -p : k)Q rp (p,t,r) 

J k = p+q 

-7 2 / d pdq qik } G^ p (q,t J r)Q p (p,t,r) 

J k=p+q 

rj pp (k .,t,r) = - f dpdqT 2 p (k,p,q)G pp {qJ J r)T 2 q {q,-p,k)Q pq (p,t,r) 

Jk=p+q 

- ~^Tfi2 J k + rfprfq ?«PjG ! “/(q,<,r)(5 p (p,< ! r) 

In this approximation, the incompressible pressure correlations appear in the expressions 
for damping of the compressible fields which take the form 


G{k , r) ~ exp[— p(k)r + 0(M 2 )} 


The appearance of M t in the response and time correlation functions is regularizing in the 
sense that time integrals 


f 


dr G(k, r )Q(k, r) 


a + 0(M 2 ) 


and the theory is therefore well behaved even if M t becomes large. This provides a more 
fundamental view of the role of pressure fluctuations and is suggestive for generalization 
of this theory to strong compressibility. 


VII. Conclusions 

TSDIA provides a systematic approach to the derivation of a transport model for 
weakly compressible turbulence. The occurrence of diffusivities for entropy and pressure 
and of shear and bulk viscosities acting on the mean flow are effects that could be antic- 
ipated by analogy with incompressible turbulence modeling. However, by incorporating 
the effects of wave motion into the description of compressible turbulence, we find Mach 
number-dependent modifications of the specific heat ratio and the speed of sound. TSDIA 
further predicts that terms dependent on gradients of turbulence quantities will appear in 
the mean flow equations. At this point, we have only given a preliminary account of these 
effects, but it is clear that they could not be predicted by other methods. 
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Nonlinear effects, analogous to normal stress effects in simple shear flow, are also 
accessible to TSDIA analysis. The lengthy analyses required to evaluate such effects can 
be expedited by symbolic calculations . 20 Other more straightforward extensions of this 
analysis are to compressible heat transfer modeling, gravitational effects, and effects of 
rotation on compressible turbulence. DIA and TSDIA provide a unified framework in 
which to analyze such coupled effects in a general and systematic manner. 

Open issues include the accurate evaluation of turbulent gradient terms and the con- 
sideration of the effects of compressibility on the two equation model. It is evident that 
these effects will be of order M 2 , but they should be evaluated for potential extension of 
this theory to strong compressibility. 

APPENDIX I. Outline of the DIA theory of passive fields coupled to incom- 
pressible turbulence 

Let a general system of fluctuating fields q i be coupled to incompressible turbulence 
by equations of the form 

q'(k,t) + A i k q k {k,t)+ f dpdq {B^(k, p, q)u~(p, f ) 

J k=p+q 

+ C i k (k,p,q)p so {p,t)}q k (q,t) + =0 (101) 

where A represents linear couplings among the q* , B and C couple the q to incompressible 
turbulence, and A couples q 1 to external fields. DIA closes all moments of the q l in terms 
of the response matrix 5 

G) = < Sq'/Bf 1 > 

where p is a small perturbation added to the right side of the equation for q J in Eq. (101). 
DIA also provides a set of evolution equations for this matrix of the general form 

<jr^-(k, t,$) + (k, f , <s) + A^kG j (k, i, s) 

+ f dr ,'(k,t,r)G*(k,r,s) = «j«(t - ») (102) 

Jo 

in which the damping factors q are 

7 ,;(k,M) = f dpdq {^(k.p,q)B' t/ (q,-p,k)Gf(p,t,5)Q M ^(q,t,5) 

J k=p+q 

+ d( k, P, q)C' (q, -p, k)Gf (p, t,s)Q p { q, t, a)} 
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Note that the damping factors 77 depend on the various couplings through their dependence 
on the response functions. 

Turbulence of the fields q l generated by some external agency is described by the 
inhomogeneous form of Eq. (101), 

9*(M) + -A^*(M) + f dpdq {B^(k,p,q) u ^°(p,t) 9 fc (q,#) 

J k=p+q 

+ CfcfoPiqK^P,*)?*^) + -4U*(M) = F'(k,f) (103) 

which has the formal solution 

<?’(M)=/ dsG)(k,t,s)Fi(k,s) (104) 

Jo 

It must be understood that Eq. (104) may only be applied to form correlations of the 
fluctuating fields q l ; DIA thus provides a statistical solution to Eq. (103), not a strict 
pointwise solution. It can be noted that DIA is the first in a sequence of successively more 
elaborate statistical approximations. 29 

Even the simplest equations of the general type of Eq. ( 102 ), such as the DIA equations 
for shear flow, 19 seem to defy analysis. A simplification which makes some progress possible 
is to treat the external field couplings perturbatively. Thus, rewrite Eq. (101) as 

q' + A)qi + B^nfq k + C>°V = -Atf (105) 

and treat the right side as a small forcing term. If the response matrix G l 0} for the left 
side of Eq. (105) can be evaluated, then it is possible to construct a formal perturbation 
series in powers of A which satisfies Eq. (105): write 

= Qo + <Zi + " • ( 106 ) 

where q' p is of order p in A , then substitute in Eq. (105). The solution of the inhomogeneous 
system Eq. (103) is Eq. (106)with 

q l 0 (k,t)= [ dsG' oj (k,t,s)F}(k,s) 

Jo 

q\(k,t) = f ds G t 0j (k,t,s)q J 0 (k,s) 

Jo 

9n(M) = / ds Go j (k,t,s)q J n-i{k,s) 

Jo 
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A particular perturbation expansion of this form was proposed by Leslie 19 to solve the 
DIA shear flow equations. This expansion is essentially a Neumann series for the response 
matrix in powers of A. TSDIA can also be understood in part as an expansion of this 
type. 

It may prove difficult to evaluate the response matrix even from the left side of Eq. 
(105). The difficulties are caused more by coupling terms like B and C than by A, which 
can often be diagonalized. A simplification which carries the procedure leading from Eq. 
(101) to Eq. (105) one step further is to treat all couplings among the q x as perturbations. 
This defines the diagonal approximation in which the response functions are evaluated from 
the simplified system 

q l + A\q l + Bi’ufq' + C'p^q 1 =0 no sum on i 

As before, the couplings can be reintroduced perturbatively; however, it will not be nec- 
essary to write the expressions generated by this expansion explicitly. Note that any 
combination of coupling terms can be treated perturbatively in this fashion. 

Applied to the system Eq. (105), the diagonal approximation will conceal an inter- 
esting physical property of compressible turbulence: the occurrence of wave motions. In 
order to reveal these motions, it is necessary to find an approximation scheme intermedi- 
ate in complexity between the diagonal approximation and the full system of governing 
equations. The approximation represented by Eq. (28) and the approximation outlined in 
Sect. VI both have this property. 

APPENDIX II. List of correlations 

A. < u°°p > 


- < u°°p > = J dk J dr 

(107) 

dp°° 

{ +G pp (k,t, r )<u?(-k,t ) 3 T (k,r)> 

(a) 

dp 00 

+ G pp {kJ,r) < uf°( k,t)U p q X (k,r) > 

(6) 

+ G pp (k,t,r) < u°°(—k,t)Upik p p 00 (k^r) > 

(c) 

dP 

+ <? pp (M,r) < ur(-M)uf(k,r) > 

(d) 


33 



0 ) 

(/) 

(g) 

(h) 
(*) 

(i) 
(*) 
(0 

Several of these terms either vanish identically, or vanish to the present order of 
approximation. As the same arguments will apply to all subsequent correlations, we give 
the argument for each term. 

Note first that G pu should be replaced by the scalar G p(t> using Eq. (31). The propor- 
tionality of G pu to k j causes several terms to vanish. 

(a), (b) both vanish to this order: correlation < u 00 ^ 00 > is zero in homogeneous turbu- 
lence; a nonvanishing contribution will require one additional gradient in A', T leading to 
a higher order term. This argument also applies to terms (e),(j),(k) 

(c) vanishes identically on integration over k because it changes sign under k — ► — k. This 
argument also shows that ( l ) vanishes. 

(d) is a gradient transport contribution analyzed below. 

(/) is a nonvanishing turbulence transport term discussed below. Other terms of this type 
are (h) and (?’). 

(g) vanishes since it is changes sign under the simultaneous transformation of wavevectors 
p — ► — p, q — ► — q, k — > — k. This transformation is admissible because it preserves the 
triangle condition k = p + q. 


+ G pp (k, t, r >7 < u~(-M)p°°(k,r) > ^ 

+ G"’(M,r) < u?°(-k,f)^(k,r) > 

+ G pp (k.t,r) f dpdq <u° 0 (—k,t)u < j°(p,r)iqjp° 0 (q,r)> 

J k=p-f q 

+ G pp (k, t, r) f dpdq < u°°(-k,t)u°°(p,r)^ F -(q,r) > 

J k=p+q OJij 

-<?J“(k,<,r)-^2 _ + d P dc l <Ui°(-k, p )p 0C (p,r)iq j p 00 (q : r) > 

- ^ dpdq < ttf , (-k,r)p 0O (p,r)^-(q,r) > 

~ G T( k ^ r ^ < «r(-M)p~(k,r) > 

+ Gf(k,f,r)<«n-k,i)C(k,r) > 
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The nonvanishing terms are (d), (/), (ft), ( 0 - Of these terms, ( d ) contains a mean flow 
gradient, and ( f),(h),(i) will depend on gradients of turbulence quantities. Note that 
although (/) vanishes to lowest order, it is possible that it makes a non vanishing contribu- 
tion at higher order in TSDIA. The resolution of this issue is left to further research. The 
terms (h),(i) are triangle integral terms of a type also discussed by Yoshizawa. 1 


B. < utf 00 > 


— < UiP > 


- / c 


dr 

dp c 


{ +C(k,<,r) <p°°(-M)-^r(Ic,r) > 

dp 00 

+ G“ p (k,tr) <p°°(-k,t)^^-(k,r) > 

+ G“ p (k. t,r) < p°°(— k,<)!7jtA; J p 00 (k,r) > 

dP 

+ G“ p (k,f,r) <p°°(- M)uf(k,r) > 

+ G“ p (k, t, r)j < p°°(— k, <)p°°(k, r) > ^ 

5p°° 

+ G“ p (M,r) < p“(-k,*)-| r (k,r) > 

+ G“ p (k,t,r) / dpdq <P°°(-k,t)«r(p,r)^ r (q,r) > 

+ G“ p (M,r) f dpdq < p^(-M)u~(p,r)i 9 j p°°(q,r) > 

Jk = p-f-q 

-G"“(k,Gr)^- / dpdq < p^-M^Cp, r^p^q, r) > 
J c 2 /? 2 J k=p+q 

1 /“ dp°° 

-G?/(M,r)-2-2 / dpdq <p 00 (-k,f)p 00 (p,r)— (q,r) > 

C ■*£ dk=p+q ^ 3 

1 dP 

- G“ u (k,f.r)^ < p 00 (-k,t)p 00 (k,r) > 


c 2 R 2 


dXj 

+ G*/(M,r) < p~(-k,t)u~(k,r) > (x^-) } 


( 108 ) 

(«) 

(ft) 

( C ) 

(d) 

(e) 

(/) 

(<?) 

(ft) 

(*) 

O') 

(ft) 

(0 


Arguments like those used in case A show that the nonvanishing terms are (c), ( h ), (j), (ft). 
Of these, only (k) depends explicitly on the mean flow. 


C. < u Vp°° > 

It will be convenient to separate the terms containing ik{ and the terms containing 
d/dXi. which will be denoted by the subscripts x and X respectively. This gives for the 
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terms containing ik{. 


- < u- V lP °° > = J dk J o dr 

{ + G“ p (M,r) < -ik iP °°(-k,t)^(k,r) > 

dn°° 

+ G“ p (k,t,r) < k,r) > 

+ G“ p (k,f,r) < -ikip^i-kJWjikjp^i^r) > 

f)P 

+ G“ p (k, t, r) < -*Ar<p°°(-k,*)ti*(k,r) > — — 

oXj 

+ G," p ( k,f,r) 7 < -ii«.“(-k,0p“(k,r) > 

< 9 d °° 

+ G“ p (M,r) < -tfc i p oo (-k,0~-(k,r) > 

+ G r “ p (M,r) / dpdq < -^ i p 00 (-k,#)u°°(p,r)^ r -(q,r) > 

-/k=p+q oAj 

+ G“ p (k,t,r) / dpdq < -^ip 00 (-k,«)w°°(p,r>' 9> p 00 (q,r) > 

«/ k= p-f~q 

- G“/(M,r)— — / dpdq < -*A; t p 00 (-k,<)p 00 (p,r) 2 9j p 00 (q,r) > 

c ^ Jk=p+q 

-Gf/(k,«,r)Ap jT dpdq < -it jp “(-k,()p~(p,r)^-(q,r) > 
-e?(k,l.r)A < -.t,p~(-k,f)p“(k,r) > 

+ G?/(k,f,r) < - t i,p”(-k,()«r( k ' r ) > (fp-) 0 } 


(109) 

(o) 

(6) 

(c) 

(<0 

( e ) 

(/) 

(?) 

(A) 

(*) 

O') 

(*) 

(0 


The nonvanishing correlations are (a), (6), (e), (/), (g), (j). Of these, only (e) depends 
explicitly on the mean flow. 

The terms corresponding to those of Eq. (109) in which iki is replaced by d/dX, are 


< u • Xxp°° > 


-hi 


dr 


{ + GHk.tr) <|£(-k,<)f£(k,r)> 


+ G“ , (k ,(,r) < ^H-k,f)l7 J ^-(k,r) > 
dry 00 

+ G“ p (k, t, r ) < {—k,t)Ujikjp°°(k,r) > 


( 110 ) 

(a) 

(b) 

(c) 
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<9n°° dP 

+ G“ p (k,f,r) < OL_(-M)u°°(k,r) > ^r- 

5e»°° 517, 

+ G“ , '(k.i.r ) 7 < ^-(-k,i)p~(k.r) > gj- 

<9»°° dv°° 

+ crik, t,r) < ^-(- k - f )-ir< k ’ r ) > 

r r^r ) 00 3 p°° 

+ G“ p (k,t,r) / dpdq < -^(-M)uf(p,r)— ^r) > 

J k=p+q aAt J 

f du°° 

+ G up {kJ,r) / dpdq < -^r(-k,t)uf(p,r)tgjP°°(q,r) > 

7k=p+q OA * 

-G“ u (kJ.r)^ / dpdq < ^ r (-k,t)p 00 (p,r)* gj p 00 (q,r) > 

J C ^ A=p+q C ' Al 


-G“”(M,r) 

-G?/(k,t,r) 


dp 

1 5 p°° . - . ~ . dP 


c 2 R 2 


/ dpdq 

7k=p+q 


p+q 

< ^(-M)j>“(k,r)> 


< if- ( ' k, ‘ )p “ <p ’ r) ^ 7 <q ' r) > 


c 2 f? 2 5A 
ap°° 


dXy 


+ G" u (k,t,r)< ^-(-k,t)«~(k,r)> (|^) C ) 


(<9 
(< e) 
(/) 
(?) 
(k) 
(*) 
O') 
(*) 
(0 


In determining which terms are to be retained, it must be noted that, unlike previous cor- 
relations, there is no leading derivative with respect to X. Therefore, terms containing up 
to two derivatives with respect to A* and T must be retained. Therefore, the nonvanishing 

terms are [c),(h),(j),(k). 


D. < u°°u > 


dr 


- < u^Uj >= J dk jf 

-G““(M,r)-^ f dpdq < uf > (—k,t)p°°(p.,r)iq m p 00 (q,r) > 

C K Jk=p+q 

i r dv°° 

-G“(M,r)^ pdq <«r(- k .‘)p“(P.’')g 3 ^(q.’-)> 

1 dP 

-Gj;(k,<,r)^<un-k,t)p~(k,r)> 

dU c 

+ GJ*(k,f,r) <«f°(-k,t)«~(k,r) > 

dp™ 

+ G" p (M,r) < uf°(-k,t)-^r(k,r) > 

<9p°° 

+ G“ p (k.t,r) < ur(-k,t)U P 2j7-(k,r) > 
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(111) 

(а) 

( б ) 

(c) 

(d) 

(e) 
(/) 



+ G7’(M,r) < uf°(— k, t)U p ik p p°° (k, r) > 

+ G“ p (k,t,r) < u°°(-M>C(k,r) > 

+ G* p (k,t,rh < un-k,0/>°°(k,r) > 

Uu\p 

+ G“'’(k.(,r)<«“(-k,()^(k,r)> 

+ <?“ p (k, t, r) f dpdq < uf°(~ k, (P, r)^ m p°°(q, r) > 

^k=p+q 

+ G“ p (k,t,r) f dpdq < u °°(-k,f)u~(p,r)|^ r -(q,r) >} 

k — p+q OJi-xfi 


(9) 

( h ) 

(*) 

(i) 
(*) 

(0 


The nonvanishing terms are (b), (d), (g), (k). Of these, only (d) depends explicitly on the 
mean flow. 


E. < u • Vu°° > 

Again writing the terms in ikj and d / d\j separately, denoting them respectively by 
subscripts x and A', the terms in ikj are 


- <(u • VjU 00 ); >= J dk J dr 

{-G™(k.t,r)-±p £_^dpdq 

-G““(k.t,r)— ^ / dpdq 
C J k=D+Q 


( 112 ) 

< -ik J u'?(-k,t)p°°(p,r)iq m p°°(q, r ) > (a) 

< —ikjU°°(—k, t)p°°(p, r)—~ (q, r ) > ( b ) 


1 Pp 

- G T^-t.r)^ < -.V“(-M)p“(k,r) > JL. 

+ G“(k.(,r) < - ! * J «r(-k,<)«r( k .’') > (l^rf 

+ G“ p (k,t,r) < -ik j uf :> (-kj)~^ r (k ,r) > 

+ G“ p (k, t, r) < -ikj U r(-k,t)U p ^-(k,r) > 

OJip 

+ G] p (k,t,r ) < -it J «f(-k,<)C7 p ik,p~(k,r) > 

+ Gf(k,(,r) < -.•i i «f(-k,t)«S(k,r) > 

C/A m 

+ (?“ p (k,t,r )7 < -?:t J uf > (-k,<)p 00 (k,r) > 

+ G] p (k,t,r) < -ik jU ?(-k,t)^(k,r) > 


( c ) 

(d) 

0) 

(/) 

( 9 ) 

( h ) 

(*') 

(i) 
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+ G“ p (k, t, r) 

f dpdq < - 2 ^uf 5 (-k,t)u^(p,r)^ m p 00 (q,r) > 

(*) 

J k=p+q 



+ G“ p (k, t, r ) 

r s dp°° 

/ dpdq < -ikjU°°{-k,t)u™(p,r) dx ( 

J k=p+q m 

q,^) >} 

(0 


The nonvanishing terms are (a), (c), (/),(/>),(*), (j), (/)• 0f these ’ ( h ) Spends explicitly 
on the mean field. 

The terms dependent on dj dXj are 


- <(u • V x u*)„ >= “ / dk J dr 

1 f du°° 

{-G““(M,r)^ + dpdq < j£-{-Kt)P°°(P,r)iqmP { ~{q,r) > 

-G>m(Mi r ) c 2^2 < ( k ’ r ) > 

+ G““(M,r) < ^-(-k,f)«~(k,r) > (|^) C 
+ G" p (M,r) < ^-(- k ^)^( k ’ r ) > 

+ G>'(k,«,r)<g(- k .wg (k ,r)> 

^7/°° 

+ G“ p (M,r) < |j-(-k,«)0 ) ,it,p~(k,r) > 

aA ; 

+ GJ'fM.r) < gj-(-k,*)u“!k,r) > 

+ G“ p (k,f,r) 7 < |£(-M)p°°(k,r) > 

< 97 , 9 ° 5 d°° 

+ GJ'fM.r) < > 

+ G" p (k,<,r) / dpdq < ^-(-k,0«m(P 5 r >'9mP°°(q,r) > 

J J k=p+q 

c f)n 00 dv°° 

+ G“ p (M,r) / dpdq < ^j-(-k,i)«~(p,r)^-(q,r) >} 

7/k=p+q J 777 


( 113 ) 

(а) 

(б) 

(c) 

(d) 

(e) 

(/) 

(*) 

(h) 

(*) 

0) 

(*) 

(0 


Again note the action of d/<9A'y on u?°(-k) onl y- The nonvanishing terms are 
(b),(d),(g),(k). Of these, (d) depends explicitly on the mean field. 
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F. < pp°° > 


- < PP°° > 


/*/ 


dr 


dp c 


{ +G"(k.(,r) <p”(-k,f)-^r(k,r) > 

+ G"(k,t,r) <p“(-k,<)^^(k,r) > 

+ G"(M,r) < p°°(— k, t)Ujikjp°°(k, r) > 

+ G pp (k.t,r) < p~(-k,«)Bf (k,r) > |£- 

+ G"(M,r ) 7 <p~(-k,i)p“(k,r) > 

+ G"(k,(,r)<p“(-k,()^(k,r) > 

+ G”(k,i,r) f dpdq <p“(-k,t) u f(p,r)^(q,r)> 
A=p+q OLAj 

+ G pp (kj,r) f dpdq < p°°(-k,#)«°°(p, r>^p°°(q, r) > 

k=p-J-q 

- Gj*(k,<,r)^j + dpdq < p°°(-kJ)p 00 (p,r)iq j p°°(q,r) > 

-G pu {k cfpdq <p°°(-M)p~(p,r)|^(q,r) > 
-Gf(k,i,r)^ <p-(-k,<)p°°(k,r) > 

+ Gf (M,r) <p~(-M) W ~(k,r) > (|^-) C } 


(114) 

(«) 

(6) 

(c) 

(<0 

(e) 

(/) 

( 5 ) 

(*) 

(*) 

0) 

(*) 

(0 


The nonvanishing terms are (a), (6), (e), (/), (p), (j). Of these, (e) depends explicitly on 
the mean fields. 


Appendix III. Approximate evaluation of turbulent correlation terms 

The evaluation of correlations which depend on gradients of turbulence quantities in 
Eqs. (108)-(112) is more elaborate. In general, it is necessary to replace p°° by u°° using 
the Poisson equation and apply the quasinormal hypothesis to close the resulting moments 
of the velocity field. Provided that the integrals converge, it is possible to determine in 
advance at least the structure of the resulting single point moments. Since it will be useful 
to have a preliminary indication of the type of model which results, the form of these terms 
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will be derived here, limiting attention to those terms -which can be evaluated at the lowest 
order in TSDIA. Shimomura’s result Eq. (52) is obtained by such higher order analysis, 
and reference can be made to Ref. 20 for the required procedure. 

A. Terms added to mean -pressure equation 
Moments from Eq. (107): 

(/) vanishes to lowest order since there are no isotropic functions of u°°. 

(h) ~ G( V. A' 4 , £ -1 , c -2 ) where the notation indicates only the total power of each factor 

in the result . 

( i ) vanishes to lowest order since it is of odd order in u°°. 

Moments from Eq. (108): 

(c) vanishes to lowest order since there is no isotropic invariant of U. 

(h) vanishes like Eq. (107/). 

(j) ~0(V,c-\K\e~ l ) 

Moments from Eq. (109): 

(a), (b) ~ 0(D/DT\ A' 5 / 2 , c -2 ) 

(/) vanishes to lowest order 

(g) ~0(V,c- 2 ,A’ 4 ,£ _1 ) 

(i) ~ Kec~ 2 

Moments from Eq. (110): 

(c) vanishes to lowest order since there is no isotropic invariant of U. 

(/?.) vanishes to lowest order 

(j) ~<9(V 2 ,A' 4 ,e-\c- 2 ) 

Term (k) is of a special type since it contains a mean velocity gradient, but its coeffi- 
cient depends on gradients of turbulence quantities. The structure of this term is 

(k) ~ [C>(V,A 3 ,£-\c- 2 )]VP. 

The mean pressure equation, with these effects added, has the form 

— + U • VP + 7 PV • U = V • (u pp VP) - (7 - l)nv • u 

dt 

c,c- 2 /v £ + c 2 o(-Ao + C 3 0(V, V, a- 4 £-V- 2 ) 
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where the constants Ci,C 2 ,C 3 can be evaluated theoretically by evaluating the relevant 
integrals. The term containing a convective derivative occurs in the higher order TSDIA 
analysis of incompressible turbulence. 

B. Terms added to mean velocity equation 

Moments from Eq. (Ill): 

( b ) vanishes to lowest order 
(g) vanishes to lowest order 
( k ) ~ K 2 c~ 2 

Moments from Eq. (112): 

(a) vanishes to lowest order. 

(e) vanishes to lowest order. 

(/) vanishes to lowest order. 

(i) vanishes to lowest order. 

(j ) ~ K 1 / 2 ec~ 2 

(/) ~ 0{V ,K 2 ,c~ 2 ) 

Moments from Eq. (113): 

( b ) vanishes to lowest order. 

( g ) vanishes to lowest order. 

(k) ~ C>(V,A' 2 ,c- 2 ) 

Term (d), like ( 1 10 A- ) depends on mean velocity gradients. It has the form 

(d) -[^(V,/! 2 ^- 1 )^ -U 

Moments from Eq. (114): 

(a),(b) ~ D/DT(K 4 £~ 1 c~ 2 ) 

(/) vanishes to lowest order. 

(g) vanishes to lowest order. 

O') ~0(V,K 9 ' 2 ,e-\c- 2 ) 

Note that this group of terms appears with an additional factor of c~ 2 and are therefore 
of higher order in M t . 
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These terms modify the mean velocity equation as follows: 

^ + u ■ vu + = V • {-2/a + *[vu + (VU) T ] S 

at R £ 

+ v(>/ uu v • u) + ivp 

R 

+ C 4 VA' 1/2 £C -2 + C 5 0(V 2 , A' 2 ,c“ 2 ) + C,0{VJ<\e- l )V ■ U 
As before, the constants C 4 ,C 5 ,C 6 are obtained by evaluating the relevant integrals. 
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